Spin-orbit coupling induced enhancement of superconductivity in a two-dimensional 

repulsive gas of fermions 
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We study a model of a two-dimensional repulsive Fermi gas with Rashba spin-orbit coupling aji, 
and investigate the superconducting instabihty using renormalization group approach. We find that 
in general superconductivity is enhanced as the dimensionless ratio ^maji/ Ep increases, resulting 
in unconventional superconducting states which break time reversal symmetry. 



There is a growing interest in materials whose inter- 
faces support a two-dimensional (2D) electron gas and 
display superconductivity, because of their novel, and 
potentially technologically useful properties such as elec- 
tronic transport, magnetism and interplay between struc- 
tural instabilities pj-yj. Due to the intrinsic breaking of 
the inversion symmetry, spin-orbit coupling is expected 
to play a role in determining the nature of the super- 
conducting state. For example, experimentally the en- 
hancement of transition temperature at LaAlOa/SrTiOs 
interfaces tracks the enhancement of Rashba spin-orbit 
coupling And while the mechanism of superconduc- 
tivity here is likely related to the electron-phonon mech- 
anism of the bulk materials 5j, such considerations moti- 
vate us to investigate the effect of spin-orbit coupling on 
the superconducting transition. 

For attractive interactions the question has been ad- 
dressed in Ref.[B-[ll- In contrast, here we consider a 
model of repulsive fermions moving in 2D and analyze 
the nature of the unconventional superconducting state 
in weak coupling. For a strictly parabolic dispersion in 
2D, without spin-orbit coupling, it is known that repul- 
sive interactions do not induce superconductivity to sec- 
ond order in the interaction, unlike in 3D where p-wave 
superconductivity is found at this order. In 2D one has 
to go to third order Q for the Kohn-Luttinger effects to 
appear. Our motivation is to understand the role of the 
spin-orbit coupling in this process, to determine whether 
it can enhance superconductivity, and to study the na- 
ture of the superconducting state. Since we treat the 
Rashba spin-orbit coupling non-perturbatively, we 
can analyze the relative values of the mean-field tran- 
sition temperatures Tc for an arbitrary value of the di- 
mensionless ratio Q = ^maj^/ Ep, where m is the (bare) 
fermion mass and Ep is the Fermi energy, measured from 
the Dirac point (see Fig. 1). In the strictest sense the 
transition in 2D is of Kosterlitz-Thouless type and at 
Tkt < Tc- However, since we are working in the weak 
coupling limit, the pairing energy scale is much smaller 
than the zero temperature phase stiffness energy and 
1-Tkt/Tc ^ Tc/Ep{l + e) < 1, justifying the approach 
presented here. 

Due to the spin-orbit interaction, the pair states can- 



not be chosen to be pure spin singlet or triplet, but 
appear as linear superposition thereof Nevertheless, 
since the Rashba model ([2]), as well as the short range 
repulsion (|3]), commute with the z— component of the to- 
tal angular momentum — Lz + Sz, we can label the 
pair states according to £, the eigenvalue of Jz- For small 
values of 9 we find that states with high values of relative 
angular momentum £ condense first, with £ decreasing as 
Q increases. For intermediate values of O we find broad 
regions of stability for £ — 4:, with dome-like dependence 
of Tc on 9, while in the limit of large 9, we find £ = 2. 
In weak coupling we show that all of these states sponta- 
neously break time-reversal symmetry. While we formu- 
late our calculation within more modern renormalization 
group (RG) approach, our results can be rederived dia- 
grammatically by summing the leading logarithms to all 
orders in perturbation theory, as has been done tradition- 
ally in treating Kohn-Luttinger effect Also, while 
our approach is similar to that of Ref. (see alsofisf), 
we use a single step RG instead of a two step RG, which 
we find more economical. 

Our starting point is the Hamiltonian for Fermions 
moving in 2D 
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where in momentum representation the kinetic energy 
(including spin-orbit coupling) is 
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and the short-range interaction energy term is 

Hint = E '^ki+k2,k3+k4CkiO-Ck2<T'Ck30-'Ck4<T.(3) 

ki . . .k4,cjcj' 

As usual, the components of k belong to the Born- von 
Karman set {27m/i} where n is an integer and L is 
the linear size of the system. Unlike in Ref.Q, we con- 
sider superconductivity for repulsive interactions, i.e. for 
M > 0, in the weak coupling limit uv2D << 1, where the 
density of states per spin in 2D for a/? = is V2d — ^ ■ 



FIG. 1: (Left) The dispersion relation. (Right) First order 
(tadpole) correction to self-energy. 



The kinetic energy term is diagonalized using the follow- 
ing transformation 
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Next, we rewrite the Hamiltonian in terms of these helic- 
ity eigenmodes. The partition function associated with % 
can be expressed in terms of the coherent state Feynman 
path integral over Grassman variables IJ] as 



PK(r)aA(T)]e- 
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where 



So ^ f dr y2 akA('r) ( + ekA - Mf ) akA(T), 
S^nt = f^(l,2,3,4)a*(l)a*(2)a(3)a(4), (6) 
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where the single particle energies are (see FigH]) 
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In the above expressions /3 = l/{kBT), fip is the ex- 
act chemical potential whose value depends on temper- 
ature T and interaction u, in such a way as to pre- 
serve average particle density. We adopt a shorthand 
expression for the multiple summations X^i 2 3 4(- ■ •) = 



f/(l,2,3,4) 



16L2 
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and a{j) = a-kjaj {tj) where aj = {/i, A, p\ and is an 
azimuthal angle in the momentum plane. 

We proceed by integrating out the high energy modes 
between the energy cutoff A and Vl <^ A about the two 
Fermi surfaces at T = 0. The expansion is organized 
by the powers of the dimensionless parameters uv2D and 
U,/A. At first order in the cumulant expansion, we find 
a correction to the chemical potential /i^ from the tad- 
pole diagram shown in FigHJ This correction is 5^f = 
-h'^{{p+) + where p± = / (|:;^ak±ak±- Such 




FIG. 2: (First row) 2""^ and 3'''' order corrections to the 4-pt 
scattering amplitude. (Second row) 4**^ order correction. For 
the 3'^'' and 4"* order terms, we display only the diagrams 
which contain logarithmic enhancement. 



negative interaction correction must be absorbed in the 
chemical potential counterterm, — fi'^p = iu((p+) + 
(/5_)) -I- 0{u^), which is positive, and which guarantees 
that the average particle density remains fixed. In gen- 
eral, we are not aware of any argument why interactions 
should not renormalize the areas of the individual Fermi 
surfaces, while of course maintaining their sum fixed, but 
to first order we find no such renormalization. 

Superconducting instability comes from second and 
higher order terms in cumulant expansion. We first find 
the renormalization of the general four fermion term and 
then we place the pairs on the two Fermi surfaces, which 
are the only processes with logarithmic enhancements. 
To second order in u, and in the Cooper channel, we 
have the following correction to the effective interaction 
action SSint = 

64L2 / ^^H$I^'^^*^^''^')"kMMa-k/.Wa-k'A(T)ak'A(T) 

kk' /jA 

where the sum over k, k' is restricted to a small window 
near the Fermi surfaces defined by indices ^ and A within 
the energy VL above and below /i^. We write 



F^,(k,k')-i;rr(k,k') + Kt(k,k'), 



(9) 



where the two qualitatively different contributions, aris- 
ing from the two 2'"^ order diagrams shown in Figl2] are 

V^^^X(k,k') = -8MA(A^++iV_)e-*"*''e*"*'''ln^ (10) 



Kt(k,k') 
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The density of states on the two Fermi surfaces are 
N± = V2D (l± In the second "particle-hole" 



contribution 



a,/3=± " ' ^ 

X (ae-'"^p - fie-'''"^) {l3e"''^+^-^' - Ae'"^"') , (12) 



where the Fermi occupation factor np(x) = 
l/(e(^-'^F')/^ + 1), evaluated in the hmit T ^ 0. 
After a somewhat tedious, but otherwise straightforward 
analysis we find that we can write 

n^A(k,k') = 2TOe-*"*''e'*>''A^A(e,cos(0k-<^k')I13) 

where A^a(6, cos((/)k — 0k)) is real. Note that under time 
reversal the helicity basis creation and annihilation op- 



erators transform as Ka\i± = =Fie*'^''a_k± and Ka 



t 

k± 



-k± 



respectively, where we used 0-k = 0k + ti". 
The above relation means that the Cooper channel poten- 
tial V^A(k, k') pairs time reversed states, as it shouldQ. 
Inspecting the form of the remaining terms in ^T0\\ as well 
as the combination A^^^ {Q , cos (f)) = ^Afj_\{Q,cos(f)) + 
|A^a(9, — cos(/)) appearing in pTt . shows that they are 
invariant under operations of the 2D rotation group. Ad- 
ditionally, since the remaining terms in the scattering 
amplitude are even under k — >■ — k, and independently 
under k' — > — k', they can be decomposed into sum over 
even angular momentum channels 



y;,^(k,k')=4me- 



cos(^(0k-0k')Ol4) 



£=0,2,4, 



where the dimensionless Fourier coefficients V^f^ are 
functions of Q and represent intra- and inter-band pair- 
ing amplitudes. 

In order to determine , we need to evaluate 
A^a(6,cos0) in Eq.^ from Eq.^. We shift p ^ 
p — |Q where Q — k — k', and transform from the 
polar coordinates to elliptical coordinates x e [l,oo), 
ip £ [0,27r) by substituting p|| = i|Q|xcos?A and p± = 
^IQIVx^ — 1 sinip. In the resulting expression ip appears 
only in cosV', so we can substitute y = cos^. For a = /3 
we then perform the integral over y first, which can be 
done in terms of elementary functions. Similarly, for 
a = —f3 we perform the integral over x first. Our anal- 
ysis is based on numerical integration of the remaining 
integral, which can be done quite fast to any desired ac- 
curacy. The final result for the antisymmetrized combi- 
nation A^'^^(e, cos 0) is shown in the Figl3] 

Next, we consider 3rd and 4th order terms in u which 
renormalize the Cooper channel. These terms can be rep- 
resented by diagrams shown in Fig [21 and used to derive 
the RG equations governing the flows of Cooper chan- 
nel couplings, which decouple in the angular momentum 




FIG. 3: Relative angle (fi — (pk — 0k' and 9 = \ma\/ Ef 
dependence of the interaction function Eas. (|ll) l T3|) and 
text below. (top) and A^i (bottom) start from ±^ at 
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dependence for finite 0, while A' 
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basis. For ^ ^ we find that the renormalized coupling 
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where . . . represents term of order which do not con- 
tain (large) logarithm as well as terms of higher or- 
der in u. If we define a dimensionless coupling ma- 
trix g^^l = ^u^niyjN^Nx V^^ and take the logarithmic 
derivative of the right hand side in ([T5]), then to, and 
including, ©(u"*), we find 



Q = ± 



(16) 



As usual, we have replaced the bare couplings by renor- 
malized couplings to the order ive are working. For 7^ 0, 
the initial condition for the above (matrix) differential 



equation is g'^^l\n=A — ^^^"^•v/A'PVaV^a . This equa- 
tion can be readily integrated by transforming into the 
orthonormalized basis for (7''[^^_^(i7) with eigenvalues 
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FIG. 4: The effective coupling appearing in the expression for 
Tc « yle-^/i^-z/i as a function of 9 = ^ma%^/EF. V2d ^ ^■ 
The dashed line at 0.0187 is the — > oo asymptote. 



where the initial eigenvalues of g^^^l\n=A, for £ ^ 0, are 
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If < for some £ or 0, then the associated renormal- 
ized coupling p7)) diverges at a scale 



(19) 



where g\^i _|_ = 2 f;^. . While the assignment between Tc 
and n* cannot reliably determine the prefactor of the ex- 
ponential term, the relative dependence on is in the 
exponential factor, which we can determine. This allows 
us to compare the dependence of the ratio of (mean- 
field) transition temperatures on aj^. For ^ = the 
equation ()16p holds as well, provided that we modify 



the initial condition to g'^^^yJ^^\vi=A = jl^^\P%W\ + 

^u^TOy^]VpVAV^|^~"\ and use the eigenvalues of this 
matrix in the Eg. dTTl) . 

To within our numerical accuracy, we find that 



V_ 



{i=Q) 



\, while K 



0, for any Q. In addi- 



tion, for <d > C(O.Ol) most dominant angle dependence is 
in while there is only very weak angle dependence 



in V-i < 0. To 0{u), .9^"°' > 0, meaning no pair- 
ing instability, and 5^"°'' — 0. To 0{u^) we find that 
^(^-0) ^ g g^j^y O > 0, due to increase in both V_[^'''' 

and V^~^^ , latter of which becomes less negative. This 
means that superconductivity resides predominantly on 
the large Fermi surface and is determined by some V^^ 
turning negative (meaning we select — in Eg. dlSp ). In 
FigUwe show the 8 dependence of the couplings for the 



9 



channel which has the highest Tc- At small value of 



0, £ is very high (see inset of Fig H]). For the intermediate 
values of 9, starting with ^ 0.005, we find the sequence 
£ = 6,4:, 6, 2, the last value of which continues to -> oo. 

Finally, we need to determine which linear combina- 
tion of the two possible zt£ states has the lowest (most 
negative) condensation energy as we go below Tc- Adopt- 
ing the arguments of Anderson and Morel [l^, we study 
this problem below Tf. within mean-field. We replace 
the full angular dependence of the pairing potential with 
just its projection on the most dominant i channel, an 
approximation which we expect to hold away from the 
boundaries separating ground states with different an- 
gular momentum. The self-consistent mean-field equa- 
tions are then solved near Tc and at T = 0. We find 
either a solution which breaks time reversal symmetry 
and fully gaps the Fermi surface(s), i.e. only one of the 
two ±£ pairing components is finite, or a solution with 
equal admixture of ±£ and with gap nodes. Comparing 
their condensation energies we find that the time reversal 
breaking solution is lower by a factor of 1.5 just below 
Tc and by e/2 w 1.36 at T = 0. For values of 6 > 0.005, 
the gap on the larger Fermi surface is much larger than 
the gap on the smaller one due to the smallness of ratio 

of V-i /V++. For smaller value of Q the two gaps may 

be comparable. 

In summary, we have studied the superconducting in- 
stability of a 2D repulsive Fermi gas with Rashba spin- 
orbit coupling. We find that due to the polarizable 
fermion background, the repulsion turns into attrac- 
tion on the large Fermi surface but not on the small 
one, giving rise to pairing there. Additional Joseph- 
son tunneling, V^_, induces pairing on the small Fermi 
surface by (weak) proximity effect. The resulting un- 
conventional superconducting states are found to break 
time reversal symmetry. While the transition tempera- 
ture is not strictly monotonic in the dimensionless ratio 
Q = ^maj^/ Ep, the general trend is that it grows with 
increasing Q. This experimentally falsifiable feature, may 
provide means for enhancement of superconductivity in 
a larger class of 2D electron systems. 

Acknowledgements: We wish to thank Prof. L. P. 
Gor'kov for useful discussions. This work is supported in 
part by NSF CAREER award under Grant No. DMR- 
0955561. 



[1] P. Zubko, S. Gariglio, M. Gabay, P. Ghosez, and J.-M. 
Triscone, Annual Review of Condensed Matter Physics 
2, 141 (2011). 

[2] A. Ohtomo and H. Y. Hwang, Nature 427, 423 (2004). 
L. Li, C. Richter, J. Mannhart, and R. C. Ashoori, Na- 
ture Physics 7, 762 (2011). 

A. D. Caviglia, M. Gabay, S. Gariglio, N. Reyren, C. Can- 
cellieri, and J.-M. Triscone, Phys. Rev. Lett. 104, 126803 



[3] 
[4] 



(2010). 

[5] C. S. Koonce, M. L. Cohen, J. F. Schooley, W. R. Hosier, 
and E. R. PfeifFer, Phys. Rev. 163, 380 (1967). 

[6] V. M. Edelshtein, Sov. Phys. JETP 68, 1244 (1989). 

[7] L. P. Gor'kov and E. I. Rashba, Phys. Rev. Lett. 87, 
037004 (2001). 

[8] V. P. Mineev and M. Sigrist, ArXiv e- prints (2009), 
0904.2962. 

[9] A. V. Chubukov, Phys. Rev. B 48, 1097 (1993). 
[10] L. P. Gor'kov and T. K. MeUk-Barkhudarov, Sov. Phys. 
JETP 13, 1018 (1961). 



[11] M. A. Baranov, A. V. Chubukov, and M. Y. Kagan, In- 
ternational Journal of Modern Physics B 6, 2471 (1992). 

[12] S. Raghu, S. A. Kivelson, and D. J. Scalapino, Phys. Rev. 
B 81, 224505 (2010). 

[13] S. Raghu and S. A. Kivelson, Phys. Rev. B 83, 094518 
(2011). 

[14] J. W. Negele and H. Orland, Quantum Many-particle 
Systems (Westvicw Press, 1998), ISBN 0738200522. 

[15] P. W. Anderson and P. Morel, Phys. Rev. 123, 1911 
(1961). 



